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ABSTRACT 

The flat-spectrum radio quasar 3C454.3 underwent an extraordinary 5-day 7-ray 
outburst in November 2010 when the daily flux measured with the Fermi Large Area 
Telescope (LAT) at photon energies E > 100 MeV reached (66 it 2) x 10~^ph cm~^s~^. 
This is a factor of 3 higher than its previous maximum flux recorded in December 2009 
and > 5 times brighter than the Vela pulsar, which is normally the brightest source in 
the 7-ray sky. The 3-hr peak flux was (85ib5) x 10~^ph cm~^s~^, corresponding to an 
apparent isotropic luminosity of (2.1ib0.2) x 10^" erg s~^, the highest ever recorded for a 
blazar. In this paper, we investigate the features of this exceptional event in the 7-ray 
band of the Fermi-L AT. In contrast to previous flares of the same source observed with 
the Fermi-LAT, clear spectral changes are observed during the flare. 

Subject headings: Galaxies: active 



Introduction 



The radio source 3C454.3, a well-known flat-spectrum radio quasar (FSRQ) at redshift z = 
0.859, has shown remarkably high activity since 2000. It has been particularly bright in the 7-ray 
band covered by AGILE and the Fermi-LAT, reaching a dail y record flux level F\E > 100 MeV] 
(Finn in units of lO^^ph cm^^s^^) of 22ibl in December 2009 (IStriani et al.ll20ld : lAckermann et al 



2OIOI ). This high flux allowed detailed analysis to be performed, making it th e best-studied blazar in 
the GeV band. Gamma-ray variability on timescales as short as a few hours (ITavecchio et al.ll2010l ) 
and a flux-doubling tim e scale of less than 3 hours for a short subflare on 2009 Dec 5 (MJD55170, 
Ackermann et al.ll20ld ) have been reported. In the LAT energy band, 3C454.3 exhibits a spectrum 
with a clear depa rture from a power-law distribution, characterized by a break ar ound 2 GeV 
(lAbdo et al.ll2009l ). This is found to be a common feature of bright 7-ray FSRQs (jAbdo et al. 
2010dl ). Little variati on of the break energy and spectra for large differences in flux states was 



observed for 3C454.3 (lAckermann et al.ll2010l ). From 77-op acity constraints, a min imum Doppler 
factor 6min ~ 13 was derived from the flux variability time (jAckermann et al.ll2010l ). and highest- 
energy photon meas urements, in a.ccord with independent measurements of 6 from superluminal 
motion observations (jJorstad et al.ll2005l ). 



In November 2010, the source displayed sustained activity at a flux of Fmn ~ 10 for several 
days before showing a fast rise to record levels of Fmn = 55, then rising to Fmn ~ 80 (as measured 
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over 6 hr-long periods). In this paper, the intraday variabihty and the associated spectral changes 
in the 7-ray band of 3C454.3 are studied and comparisons are made with the findings obtained 
from earher major flares. In Section 2, observations and analysis of Fermi-LAT data from 2010 
September 1 to December 13 are presented. Results are presented in Section 3 and discussion is 
given in Section 4. A flat ACDM cosmology with Hq = 71km s~^ Mpc"\ = 0.27 and f2A=0.73 
is used in this paper. 



2. Observations and analysis 



The analysis performed for this paper is very similar to that reported in lAckermann et al 



(|2010l ). to which we refer for details. The data presented in this paper are restricted to the 
100 MeV-200 GeV range and were collected from MJD55440 (2010 September 1) to MJD55543 
(2010 December 13) in survey mode. 

Spectral analyses were performed by fitting the spectra with multiple different models over the 
whole energy range covered by the LAT at E > lOOMeV. The spectral forms considered are a broken 
power law (BPL, N(E) = Nq{E / E^^^^^Y^^ with i = I \i E < E^^^^j^ and i = 2 li E > E^^^^j^, 
a log-parabola function {N{E) = Nq {E/Epy^-I^ \og{E/Ep)^ where Ep is fixed at 1 GeV), a power 
law with exponential cutoff function (PLEC, N{E) = Nq {E / Ep)~^ exjp{— E / E ^^^^jj)) , and a PL 
model over equally spaced logarithmic energy bins with T kept constant and equal to the value 
fitted over the whole range. 

Source variability was investigated by producing light curves with various time binnings (3hours, 
6 hours, 1 day, 1 week) and over different energy ranges (E>100 MeV, E>1 GeV, E=0.1-l GeV). 
Although the actual spectral shape exhibits definite curvature, light curves were produced by mod- 
eling the spectra in each time bin as a simple power law (PL) over the considered energy range, 
since the statistical uncertainties on the power-law indices are smaller than those obtained from 
BPL fits. In order to minimize spurious correlations between integrated flux and F, the fluxes 
Fe> Eq were also computed above the "decorrelation energy" £"0 where thi s correlation is minima l. 



For the 2009 December and 2010 April flares, Eq was found to be 163MeV (lAckermann et al.ll2010l ). 
The same value has been used here for consistency. The estimated systematic uncertainty on the 
flux is 10% at 100 MeV, 5% at 500 MeV and 20% at 10 GeV. The energy resolution is better than 
10% over the range of measured E^^^g^^. 



3. Results 

Figure [1] (top panel) showing the historical Fioo light curve illustrates the spectac- 
ular rise in activity of 3C454.3 over the years. It is evident that the December 2010 
outburst dwarfs any previously recorded flares. The second panel displays the light 
curves with time binnings of 1-day (open circles) and 1-week (filled circles) during the 
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outburst period. A 13-day-long plateau precedes the 5-day-long flare, confirming the trend pre- 
viously observed in the July 2008 and December 2009 flares, but it is longer in duration and higher 



in flux than those in previous flares (lAckermann et al.ll2010l ). The onset of this plateau is clearly 



accompanied by a weak but significant hardening of the spectrum (F changes from 2.50 it 0.02 to 
2.32 lb 0.03), as observed on a weekly time scale in the bottom panel of Figured) The daily flux 
decreases by a factor of about 3 in 4 days at the end of the flare. The flare is followed by a slowly 
decaying activity around -Fioo = 20. Different time periods labeled pre-flare, plateau, flare and 
post-flare in Figured] are considered in the following. 

The light curves for Fiqq with 6-hr and 3-hr time binnings (focusing on the flare period) are 
shown in the upper panel of Figure [2j The F[£' > 1 GeV] light curves with 6-hr and 3-hr time 
binnings are given in the second panel of Figure EJ The corresponding evolution of F is plotted in 
the third and bottom panels. As can be seen, the major flare lasts for about 5 days. In the Fioo 
light curves, it seems to comprise three to four subflares, with the flux peaking during the last one. 
The rise time of the MJD55516.5 flare is 12 hours for a factor of 4 increase in flux (i.e., a doubling 
time of 6 hours). This is shorter than the previous fastest relative flux variation claimed in the 
GeV band for a major flare, which was from PKS 1502+ 106, when an incr ease by a factor of 3 



in 12 hours (i.e., a doubling time of 7.5 hours) was found (jAbdo et al.ll2010d ). In the upper panel 
of Figure [21 a fit consisting of a slowly varying background and four faster temporally evolving 
components was performed between MJD55516.5 and MJD55522 for both the 6-hr and 3-hr light 
curves. Each component is assumed to be fit by a function of the form 

F = 2Fo(e(*»-*)/^'- + e(*-*«)/^0"' (1) 



(jAbdo et al.ll2010bl ). where Tr and Tf are the rising and falling times, respectively, and Fq is the 
flux at to representing approximately the flare amplitude. With the rise time Tr set equal for all 
subflares, and likewise for the fall time Tj, we flnd that Tr = 4.5 it Ihr and = 15 it 2hr gives 
a good fit to both the 3-hr and 6-hr light curves. The highest-energy photon collected during the 
MJD55516-55522 period within the energy and inclination angle-dependent 95% containment angle 
around the source position has Emax = 31 it 3 GeV and was detected at MJD55521.46. Its detection 
time is depicted with an arrow in the upper panel of Figure O 

Significant differences between light curves for fiuxes in the E=0.1-l GeV (Fo.i_iGev) and 
E> 1 GeV (-FiGev) ranges are observed (see the second panel of Figure [2|). The peak of the first 
subflare occurs approximately 15 hr after Fo.i_iGeV has leveled off, demonstrating clear spectral 
variability. This behavior is confirmed by a progressive decrease of F (spectral hardening) from 
F « 2.35 to F 2.1 as the subflare develops beyond the 100 MeV peak (black points in the third 
panel of Figure [2]). Overall, the -FiCeV light curve shows sharper structures than the -Fo.i_iGeV light 
curve. A clear difference is also observed in the decaying stage, with the high-energy component 
starting to fade away later than the lower-energy component. In the lower two panels of Figure O 
the spectral indices obtained over the two restricted energy ranges are shown as well. As Fo.i-i GeV 
levels off in the first subflare while Fi GeV keeps rising, the 0.1 — 1 GeV index is fairly constant. 
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indicating that the hardening is hmited to the range above 1 GeV. This hardening is confirmed by 
the evolution of T measured in the > 1 GeV range (blue points in the bottom panel in Figure [2]). 

The pronounced spectral evolution observed during the flare can be further investigated by 
plotting r versus the flux above Eq = 163 MeV. This is done, with a 6-hr binning in Figure [3] for 
two consecutive time periods covering approximately the first and second half of the major flare 
and two different photon energy ranges: E > 0.1 GeV and E = 0.1 — 1 GeV. For orientation, the 
points associated with the earlier times have labels corresponding to those given in the third panel 
of Figure [21 In addition, 4-day averages obtained during the plateau period are displayed as blue 
squares. In contrast to the 2009 December-2010 April flares for which no clear pattern was found 
in the data, a clockwise pattern is observed for the second period. The reduced Xr ^ with a 
constant F are 37.6/8 (P = 8.9 x 10"^, ~ AAa) and 35.4/8 (P = 2.3 x 10-^ ~ 4.2cj) for the first 
and second periods, respectively (for the E > 0.1 GeV case). For the first period, a flux increase by 
a factor of 4 is accompanied by an essentially constant F (or one becoming weakly harder). This is 
followed by a clear hardening of the spectrum at constant flux, with F changing from 2.24 it 0.06 
to 2.11 ±0.04 in 12 hours. 

Such a hard lag can also be observed in the top right-hand panel of Figure O where a hardening 
by A F = 0.42 ± 0.13, associated with the decaying stage of the flare, occurs over 2.25 days. The 
spectrum softens fairly quickly afterward, with F changing from 2.0 it 0.1 to 2.31 it 0.07 in 6 hours. 
This behavior may be driven by cooling. During the flare, the electron energy distributions may 
reflect the alternative dominance of acceleration and cooling processes. 

Restricting the analysis to the E=0.1-l GeV range (bottom panels in Figure [3|) produces a 
pattern similar to the E>0.1 GeV case for the second period, but somewhat different for the first 
period, although the patterns are less clear due to larger statistical uncertainties. The rise in flux 
is accompanied by a pronounced hardening in this energy range followed by a state of essentially 
constant spectral hardness evolving into a slow softening. This behavior confirms the conclusions 
obtained in the context of Figure [2l Due to insufficient statistics, no clear pattern of the photon 
spectral index above 1 GeV vs flux can be observed with 6-hr time binning. 

Figure 13] shows the vFiy distributions of the Fermi-LAT data for the four different time periods 
delineated in Figured! These distributions have been fitted with BPL (solid), log-parabola function 
(dashed), and PLEC (dashed-dotted) functions. The parameters of the different fits can be found 
in the Tabl^H As the likelihood method does not provide an absolute goodness-of-fit measure, 
the of the u Fy data points for the different functions have been calculated. For the pre-flare 
and plateau periods, both BPL and PLEC give fits of similar quality, while the log-parabola fit 
is notably worse. The PLEC function is preferred for the post-flare period. None of the tested 
functions provides a very good fit to the energy distribution in that period, which may be a result 
of the significant spectral evolution during the flare. The v Fi, spectra obtained over time intervals 
where the four subflares alternatively dominate are consistent for the first three subflares in terms 
of curvature, while a significantly harder spectrum is observed between MJD55520. 0-55521. 5. In 
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that interval, the PLEC fit gives E^^^^^jj =8.3±1.7 GeV. A total of 10 photons with E > 10 GeV 
(out of 12 detected during the entire 5-day flare) were collected in that 1.5-day time lapse (second 
panel of Figure [2]) . The variation of Ef^^^^j^ and E^y^^^jj with flux are displayed in the inset of 
Figure HI As already found during the 2009 December and 2010 April flares, no strong evolution 
of either -E-^^g^^ or E^^^^jj is found. E^^^^j^ remains constant within a factor of ?a 2 while the flux 
varies by a factor of ~ 40. 



4. Discussion 

During its five-day outburst from 2010 November 17 to 21 (flare interval in Figure [1]), 3C454.3 
was the brightest GeV 7-ray source in the sky, with a flux Fioo = 66 it 2 on 2010 November 18-19. 
Prior to the flaring phase, the Fermi-LAT light curve displays a 13-day long flux plateau preceding 
the major outburst. The onset of the plateau is marked by a rapid (< 1) day flux increase by a 
factor of ~ 2. This feature appears to be a c haracteristic behavior indicating that 3C454.3 is about 



to flare, as noted in lAckermann et al.l (j2010l ). In the December 2009 flare, the plateau lasted for 6 
days at a level of Fioo ~ 10 before flaring to a daily flux of ~ 22, while for the April 2010 outburst, 
it lasted for 7 days at a level of Fioo ^ 7 before reaching a peak flux of ~ 16. The spectrum 
hardens slightly from the pre-flare to the plateau preceding the giant flare. Spectral hardening 
and clustering of photons with E > IQ GeV is also seen in the decaying stage of the gamma-ray 
outburst at MJD55520. 0-55521. 5, which could point to the presence of radiating hadrons or 77- 
absorption effects. In the former case, protons require additional time to accelerate and cool while 
the < 1 GeV flux, if due to rapidly cooling electrons, would decline more rapidly. In the latter case, 
the emergence of the hard component could occur if the radiating plasma becomes optically thin to 
77-absorption, either due to a larger bulk Lorentz factor or increased size of the radiating plasma. 

The features of the giant flare can be compared to those of the two earlier, fainter flares 
(December 2009 and April 2010) that have been carefully investigated in the LAT energy band 
(the different observation mode used during most of the July 2008 flare provided poorer-quality 
data). The overall light curves show similarities (presence of a preflare plateau, main flare lasting 
a few days, several week-long fading period). The rise time of the November 2010 flare is about 
half that of the December 2009 flare (4.5 hr vs 8.9 hr). Whereas the latter showed indication of 
"flickering" activity on timescales as short as 3 hours above 100 MeV, this effect is not clearly 
present here, as demonstrated by similar 3-hr and 6-hr light curves in Figured! For the first time, a 
significant temporary hardening of the spe ctrum leading to F :^2.1 has been observed for 3C454.3. 



Note that in the first LAT AGN catalog (jAbdo et al.l l2010al ). less than 2% of FSRQs are found 



with 11-month averaged F <2.1. The moderately h ard spectrum during t he large luminosity flare 



deviates from the trend seen in the blazar divide (jGhisellini et al.l |2009| ) , where the most 7-ray 
luminous blazars generally have F > 2.5. Despite the overall spectral variation, the energy cutoff 
remains essentially unchanged as observed in earlier flares (Figure U]) . Interestingly, several-day 
long spectral variations are also observed during the postEare period (beyond MJD55524 in Figure 
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[T]). No such effect was found in previous flares despite sufficient measurement statistical accuracy. 
Significant spectral hardening at the end of the main flare was not seen either. 

The minimum Doppler factor Smin can be numerically evaluated from 77-opacity constraints. 
Prom Swift-XRT public data, the total energy flux in the 2-10 keV range is 0.8 x 10^^'^ergcm~^s~^ 
and the photon-number index is 1.70. Correlated X-ray and GeV variability supports the assump- 
tion that the 7-rays are made co-spatially with the X-rays. At the time that the 31-GeV photon 
was observed, i„ar=0.4 d (taken as ln(2) x Tf in the fading phase of the flare), giving 6m.ip . = 16, 
which is somewhat larger than the value of 5min = 13 found by lAckermann et al.l (|20ld l. The 
estimated comoving size of the emission region is R' = ct^ar ^min/i^ + -2) ~ 3 x 10^^ cm. For the 
31-GeV photon, the optical dep th to pair production by the extragalactic background light (EBL) 
is ~1 for the high-EBL model of lStecker et al.l ()2006l ). so no absorption constraints are provided by 
these data. 

An upper limit on the optical depth T^'y{Emax) ~ 2 arising from the interaction of 7 rays with 
broad-line-region (BLR) photons can be obtained by comparing the flux measured at Emax=31GeV, 



with the flux extrapolated from lower energy. Assuming a BLR luminosity of 3 x 10^^ erg s ^ (jPian et al 



20051) and adopting a BLR size rRT,R ~ 6 x 10 cm ~ 0.2 pc from reverberation mapping 



( Kaspi et al. 2007 : Bonnoli et al. 201 ll ). we c alculate t^^ (z pm)-. where Zem is the distance of the 



emitting blob from the black hole (following iReimed 120071 ) . Assuming that the BLR clouds are 
distributed between 0.01 pc and vblr, the condition T^^{Emax) = 2 is satisfied for Zem — 0.14 pc. 
Therefore the emission region must have been located either close to the outer boundary of the 
BLR or beyond at the time of emission of the 31-GeV photon. 

The asymmetry of the time profiles derived for the subflares can be pro duced by accelera- 



tion/radiativ e losses or light-travel effects in quasi-spherical emission regions (jSikora et al.l 12001 



Dermerjl2004l ). Assuming a spherical geometry of the emission region, then r < 2cr^ t^ar/ (1 + -2) ~ 
0.1 pc for a bulk Lorentz factor F = 5min = 16 and a variability time scale t^ar — 0.4 d. If the jet 
opening angle 6j <C 1 /Smin — 3° , then this estimate can be compatible with the location estimated 
above. 

The 3-hr peak Fioo is 85±5, corresponding to an apparent isotropic 7-ray luminosity = 
IO^OLso erg s'^ with L50 = 2.1 ± 0.2, assuming a spectral shape for the flare as given in Table 1. 
This is ~ 4 times the luminosity of PKS 1 622—297 (L50 — 0.5 with the current cosmological model) 



during its 19 95 flare (IMattox 



for a blazar. (Foschini et al 



et a 



( 2011 



19971 ). making this the largest 7-ray luminosity ever observed 
derived a luminosity of L50 ~ 3.0 during this flare by assuming 
a power-law shape and considering a short 4.7 ks time interval where Fioo ^ 100.) During the 



5-day flare, Fioo = 43 ± 1 imp Hes L^n 



the range of (0.5-4) x 10^ M0 (jBonnoli et al. 



1.0. The black hole mass for 3C454.3 is estimated to be in 
Gu et al.l[200ll ). so LEdd ~(0.6-5) x lO'''^ erg s^^. 



2011 



In order for the time-averaged flare luminosity to be lower than LEdd) &j ^ 2°-6°, and a beaming 
factor (1 — cos6j)~^ > 200-1700 is implied. For highly efficient 7-ray production, t he absolute jet 



power is comparable to the disk luminosity of =6.75x10''^ erg s ^ estimated in iBonnoli et al 
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(|201lh . 



This flaring episode marks a record among AGNs and, indeed, all non-GRB sources for its 
rate of change in apparent luminosity, L^/At. Using a 6-hour variability time scale, then L^/At ~ 
lO^'^Lso erg s~^/(10^t4 s) ~ l O'^^L^n / td ergs~^, notin g that 6 hr/(l + 2;) ~ 10'^ s. By comparison, the 



giant flares of PKS2155-304 (lAharonian et alJ^nOTl ) reached only < lO'^'^erg s-V300s ^ 3 x 10"^^ erg 



s~^ (in the TeV regime). This value now greatly exceeds L^dd/iRs/c) ~ 1.3 x 10^^ erg s^^ given 
by the ratio of the Eddington luminosity and the light crossing time across the Schwarzschild 
radius of a black hole, and strongly violates optically thin, Eddington- limited accretion scenarios 
(|Elliot ShaDirJl974l ). showing that such a condition is unlikely to apply to the highly asymmetric 
disk/jet system of a blazar. 
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Fig. 1. — Top, main figure: 
LAT since launch. Inset: 



Daily light curve of 3C454.3 measured with the Fermi- 
Historical light curve. Black points are from E GRET 

Bot- 



2OI0I). 



(jHartman et al.l 1 1993 ) and blue points are from AGILE (jStriani et al. 

torn: Light curve of the flux Fioo (top) and F (bottom) for a 103-day period including both the 
slowly increasing flux phase, the plateau, the flare and the post-flare. The open and filled symbols 
correspond to daily and weekly averaged quantities respectively. Error bars are statistical only. 
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« 0.1 — 




55514 ■ ' 55516 ' ' 55518 ' 55520 ' 55522 ' 55524 

Time (MJD) 



Fig. 2. — Top panel: light curve of the flux above 100 MeV. Open and filled symbols correspond 
to 3-hr and 6-hr averaged quantities respectively. The lines correspond to the results of a five- 
component fit (four subflares and a slowly- varying background) as described in the text, using 
the 6-hr data. The arrow depicts the detection time of the 31-GeV photon. Second panel: flux 
Fo.i_iGeVi multiplied by a factor of 0.05 (red); flux FiceV (blue). The blue histogram represents the 
times of detection of > 10 GeV photons. Third panel: T measured at E>100 MeV (black) and 
E=100 MeV-1 GeV (red). Bottom panel: T measured at E>100 MeV (black), E=100 MeV-1 GeV 
(red) and E>1 GeV (blue). 
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Fig. 3. — r vs. flux above the decorrelation energy for two different periods of time during the 
flare, obtained with a 6-hr binning. Photons with energy above 0.1 GeV (top) or in the 0.1-1 GeV 
range (bottom) were used. In the left panels, the blue points correspond to 4-day averages in the 
"plateau" period preceding the flare. 
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Fig. 4. — Spectral energy (uF^) distributions for four different time periods (pre- flare: magenta, 
plateau: green, flare: blue, post-flare: red), along with the fitted BPL (solid), log-parabola (dashed) 
and PLEC (dashed-Sdotted) functions. The inset displays -Ebreak (red) and EcutoS (blue) as a 
function of flux for the different periods considered here (filled symbols) and for historical flares 
(open symbols). 
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period 


^100 

(10-6 ph cm-^g-i) 


Luminosity 
(lO^S erg s-i) 


Ti 




nTPn Jc 

(GeV) 


AL 


')Ct 


1 


3.52±0.08 


7.8 


2.34±0.02 


2.95±0.07 


1 o+°-i 


-31.6 


6.3/8 


2 


11.2±0.2 


26.3 


2.28±0.02 


3.00±0.10 


9 O+0.3 


-18.1 


6.5/9 


3 


43.0±0.6 


105.8 


2.15±0.01 


2.81±0.05 


1 7+0-1 


-74.5 


45.9/9 


4 


20.2±0.3 


45.5 


2.29±0.02 


3.20±0.10 


9 q+0.3 

^••-'-0.3 


-44.4 


16.6/8 


period 


^100 

(10-6 ph cm-^s-i) 


Luminosity 
(lO^S erg s-i) 


a 


/5 




AL 




1 


3.45±0.07 


7.7 


2.61±0.03 


O.llibO.Ol 




-26.3 


19.6/9 


2 


10.9±0.3 


26.6 


2.39±0.02 


0.06±0.01 




-13.8 


12.5/10 


3 


41.7±0.7 


103.5 


2.36±0.02 


O.llibO.Ol 




-73.7 


43.9/10 


4 


19.1±0.4 


44.4 


2.49±0.02 


0.12±0.01 




-37.4 


13.0/9 


period 


Fioo 

(10-6 ph cm-^s-i) 


Luminosity 
(lO^S erg s-i) 


r 




^cutoff 
(GeV) 


AL 


X? 


1 


3.5± 0.1 


7.7 


2.30 ± 0.04 




5.0±1.0 


-24.1 


9.4/9 


2 


11.1±0.4 


26.0 


2.23±0.03 




11.0±2.4 


-18.3 


7.2/10 


3 


42.8±1.0 


102.1 


2.09±0.02 




6.2±0.7 


-85.5 


22.8/10 


4 


20.0±0.5 


44.8 


2.21±0.02 




5.9±0.8 


-46.7 


6.3/9 



Table 1: Parameters of the BPL (Fi, ^break^^ log-parabola (a, /3) and power-law+exponential 
cutoff (r, ^Quiojj) functions fitted to the spectra for the different periods considered in Figur^H 
AL represents the difference of the logarithm of the likelihood with respect to a single power-law 
fit, and Xr represents the reduced chi-squared of the ly Fy data points for the different functions. 



